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FOREWORD 


This report presents the second year results of an investigation being con- 
ducted under Contract NAS8-26363 for NASA George C. Marshall Space 
Flight Center under the technical direction of the Aero-Astrodynamics 
Laboratory, Dynamics and Control Division. Dr. S. Winder is the tech- 
nical monitor. The study is being performed by Convair Aerospace 
Division of General Dynamics under the direction of Mr. R. Huntington, 
project leader. 

The author is indebted to Mr. H. Riead for his assistance in performing 
the numerical analysis. 
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SUMMARY 


A totally reusable space shuttle configuration has been analyzed during 
ascent flight to determine its response to atmospheric turbulence. Re- 
sponses in the form of booster and orbiter body accelerations and booster 
wing root shear, bending moment and torque were obtained due to random 
and quasi-square -wave discrete turbulence. The configuration was also 
analyzed with booster aerodynamic surfaces removed to simulate an 
expendable booster. Symmetric and antisymmetric analyses were per- 
formed. Propellant sloshing, gust penetration, and automatic control 
system effects were included. 

The steady-state aerodynamic method of Woodward was used to derive 
the generalized aerodynamic forces using the standard quasi-steady 
assumption. Aerodynamic interference effects between adjacent wings 
and bodies were found to be significant. Total normal force and moment 
data computed by the Woodward method agreed well with wind tunnel values. 

It was found that the symmetric responses were generally higher than the 
antisymmetric ones. The stability augmentation system tended to lower 
the booster accelerations in the symmetric case , while increasing the 
orbiter accelerations. 

Loads due to the 9 m/s quasi-square-wave gust were higher than the 3 ct 
random turbulence loads , indicating that the discrete gust design criterion 
may be conservative . The elastic portion of the response accounted for 
about 15% of the total wing load in the discrete gust analysis, while in the 
random case the elastic effect was small. 
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SECTION 1 


INTRODUCTION 


This report presents the results of the second year's work under Contract NAS8-26363, 
Aeroelastic Effects on Shuttle Vehicle Dynamics . The first year of the study is docu- 
mented in References 1 and 2. 

During the first year, a method was developed for determining the response of space 
shuttle to atmospheric turbulence. The objective of the second year's effort was to 
apply this technique to the shuttle configuration during ascent flight. 

The analysis method accounts for propellant sloshing, gimbaled engine dynamics atid 
stability augmentation system (SAS)/elastic vehicle coupling. The problem is formu- 
lated in the frequency domain to permit the use of existing unsteady or quasi-steady 
state aerodynamic theories based on harmonic motion for computing generalized aero- 
dynamic forces. Either discrete or random turbulence can be treated. For the 
random case , statistical outputs are generated relating vehicle loads and acceleratiotis 
to level of turbulence. A Fourier series is used to approximate the discrete gust 
shape. Discrete responses are in the form of time histories. 

A finding of the first year's study was that properly predicting flow interference, 
thickness, angle -of -attack, and body aerodynamics is more important for space shuttle 
than including unsteady aerodynamic lag effects for the low to moderate reduced- 
frequency range of interest in the turbulence response problem. It was recommended 
that the Woodward steady-state aerodynamic method (Reference 3) be applied in future 
space shuttle turbulence response studies using the standard quasi-steady state 
assumption which includes aerodynamic damping due to plunging velocity. 

In the study reported here, the space shuttle configuration shown in Figure 1-1 was 
analyzed at Mach numbers of 0. 8 and 1.2, which are in the region of maximum dynamic 
pressure on the ascent trajectory (Figure 1-2). Symmetric and antisymmetric 
responses to random and discrete turbulence were obtained. The turbulence models 
were as specified in Reference 4. The Woodward aerodynamic method was used to 
compute the generalized forces. To simulate a non-flyback booster , calculations were 
also made for the configuration shown in Figure 1-1 with all booster aerodynamic 
surfaces removed. Propellant sloshing effects were included. The SAS was included 
in one symmetric condition, A summary of conditions analyzed is given in Table 1-1. 

Section 2 contains detailed discussions of the modal and aerodynamic analyses , and 
results of the turbulence response calculations. Conclusions and recommendations 
for further work are presented in Section 3 . 
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Figure 1-1. Study Vehicle 
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Figure 1-2. Ascent Trajectory 
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Table 1-1. Conditions Analyzed 


Condition 

No. 

Mach 

No. 

Booster Aero 
Surfaces 

Boundary 

Condition 

SAS 

Response 

Type 

1 

1.2 

On 

Sym. 

Off 

Discrete 

2 

1.2 

On 

Sym. 

Off 

Random 

3 

1.2 

On 

Sym. 

On 

Random 

4 

1.2 

Off 

Sym. 

Off 

Random 

5 

1.2 

On 

A/S 

Off 

Random 

6 

1.2 

On 

A/S 

Off 

Discrete 

7 

1.2 

Off 

A/S 

Off 

Random 

8 

0.8 

On 

Sym. 

Off 

Random 

9 

0.8 

Off 

Sym. 

Off 

Random 

10 

0.8 

On 

A/S 

Off 

Random 

11 

0.8 

Off 

A/S 

Off 

Random 
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SECTION 2 


TECHNICAL DISCUSSION 

Significant effort was put forth in deriving the input data for the turbulence response 
analysis. Therefore, Sections 2.1 and 2.2 present details concerning the modal and 
aerodynamic analyses, respectively. A detailed discussion of the turbulence response 
analysis, along with significant results, is presented in Section 2.3. 

2.1 MODAL ANALYSIS 

The NASTRAN computer program (Reference 5) was used to obtain the structural 
vibration modes for the mated space shuttle configuration. Hie over-all structural 
idealization is shown in Figure 2-1. The orbiter body is represented as a flexible beam, 
and the orbiter aerodynamic surfaces are rigid. The booster body and aerodynamic \ 
surfaces are elastic, represented in considerable structural detail. By use of the 
Guyan reduction (Reference 5), the dynamic degrees of freedom were reduced to 107 
for the symmetric analysis and 106 for the antisymmetric case. 



Figure 2-1. Over-all Structural Idealization 
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The modal frequencies are summarized in Table 2-1. The mode shapes for t = 75 
seconds are shown in Figure 2-2 for the symmetric case. Figure 2-3 presents the 
corresponding antisymmetric modes. 


Table 2-1. Summary of Modal Frequencies (Hz) 


Type 

t = 60 sec 

t = 75 sec 

Symmetric 

1.07 

1.09 



1.17 

1.17 



1.48 

1.48 



2.05 

2.06 



2.18 

2.19 



2.49 

2.49 



2.80 

2.89 

1 


3.02 

3.04 

Antisymmetric 

0.87 

0.87 



1.26 

1.28 



1.77 

1.77 



2.37 

2.52 



3.19 

3.19 



3.45 

3.82 



4.28 

4.32 


The mode shapes at t = 60 seconds were assumed to be identical to those at 75 seconds. 
The generalized masses were corrected for the heavier mass distribution by the 
equation 


n 

%0. = %5. Am i ^i (1) 

J J Pi j 

where is the generalized mass of mode j at time t, Amj is the change in mass at 
.th ^ 

the i structural grid point, 0$. is the deflection at grid point i due to mode j, and the 
summation is over all grid points, n. 

The frequencies were corrected by the equation 
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( 2 ) 




Figure 2-2. Symmetric Mode Shapes (t = 75 sec) (Sheet 1 of 2) 




Mode 7, 2.89 Hz Mode 8, 3.04 Hz 

Figure 2-2. Symmetric Mode Shapes (t = 75 sec) (Sheet 2 of 2) 





The effect of removing the booster aerodynamic surfaces was accounted for in a 
similar fashion. The mode shapes and frequencies were assumed to remain constant, 
and the generalized masses were adjusted according to Equation 1. Modes consisting 
almost entirely of booster aerodynamic surface motion were eliminated (symmetric 
mode 6 and antisymmetric modes 1, 5 and 7). 

2.2 GENERALIZED AERODYNAMIC FORCE DETERMINATION 

The steady-state aerodynamic method of Woodward (Reference 3) was applied to the 
space shuttle configuration to generate aerodynamic forces. The justification for 
using steady-state aerodynamics instead of unsteady is given in Reference 1. T ° 
summarize briefly: the Woodward method accounts for wing-wing, wing-body and body- 
body interference, angle -of -attack, thickness and other flow effects not accounted for 
in most unsteady theories. The reduced frequency range of interest is low enough in 
the shuttle turbulence response problem that aerodynamic lag effects are relative y 
small. By using steady-state aerodynamics, the computation of aerodynamic forces 
is made once for each Mach number instead of for each value of reduced frequency as 
is required in the unsteady case . 

Assuming harmonic motion, the quasi-steady complex downwash at point j due to a unit 
amount of the k^ 1 generalized coordinate is 

= v “jk - 1 “ h jk t3> 

where the first term is the downwash due to the flight velocity, V, and the modal angle 
of attack, a; and the second term is the downwash induced by the surface plunging 
motion of amplitude h and frequency W . 

The generalized aerodynamic force consists of two parts; that due to turbulence, and 
that due to vehicle motion. Specifically, the generalized aerodynamic force acting on 

mode r is given by 

Qr = % Qrs + Tf Qrf ^ ^ 

s=l 


Response Gust 

where QL is the r th generalized aerodynamic force due to a unit value of the vehicle 
displacement in mode s , q f is the gust amplitude , and Q rf is the r* generalized aero- 
dynamic force due to a unit gust. is defined by the equation 


% S (U)=J h r (x,y,z) Ap s (x,y,z,W) dS 
S 
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where AP S is the net pressure acting on the vehicle due to oscillation in mode s , and 
the integration is over the surface of the vehicle . 


The generalized force due to the unit sinusoidal gust is given by: 
Q r f(tO) = / h^x.y.z) Apj>(x , y , z , 40 ) dS 


( 6 ) 


where Apj is the net pressure over the vehicle surface due to the downwash produced 
by a continuous sinusoidal gust wave traveling across the vehicle. This sinusoidal 
gust velocity may be expressed as: 


W, 


.icot 


'f - Qfe 

= (cos cot + i sin cot) ( 7 ) 

The downwash on an oscillating airfoil in a flow of velocity (V) is given by: 

w(t) - V a(t) -h(t) < 8 ) 

where a is the angle of attack and h is the plunging velocity. For harmonic motion, 
the complex downwash may be written: 


w = V a - i coh 
Equating (7) and (9) yields: 

W f = w 

— qf 
a = ~ cos cot 

h = - ^£-sin cot 
co 


O) 


( 10 ) 

( 11 ) 


If xj is the distance from a gust reference axis to the downwash at point i on the vehicle 
and k is the reduced frequency 


(-t) 


( 12 ) 


then Equations 9 and 10 become: 


a i = cos k xj 


( 13 ) 
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( 14 ) 


h. = - 
1 


qfb . . - 
fk smkx i 


where 



Equations 13 and 14, when substituted for a and h in Equation 3, give the downwash 
distribution over the vehicle due to unit turbulence. 

Having the uni t downwash distributions , the Woodward method is used to derive the net 
pressure distributions required to solve Equations 5 and 6. Essentially, a matrix of 
influence coefficients, [F] , is generated relating pressure to downwash in the form 

|pj=[F]{w) >> 

Woodward represents the vehicle as a system of source, doublet, and vortex singu- 
larities. The wings and bodies are divided into panels defining the boundaries of planar 
singularities. Line singularities distributed along the body centerline simulate body lift. 

The aerodynamic panel representation employed here, shown in Figure 2-4, consists of 
91 wing and 126 body panels. The body centerline points locating the line singularities 
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(not shown) were generally distributed per each ring of body panels and continued 
forward to the nose of each body at about the same spacing. 

The booster wing, canard, and tail panels were removed to represent the configuration 
without booster aerodynamic surfaces. 

There is no convenient way to display the aerodynamic influence coefficient data gen- 
erated in this study. There were some interesting results, however, which will be 
discussed briefly. 

Figure 2-5 shows the nondimensional aerodynamic center location and the total normal 
force coefficient derivative as a function of Mach number for the booster in the presence 
of the orbiter. The continuous curves are based on wind tunnel data as reported in 
Reference 7. The points represent the Woodward results at Mach number 0.8 and 1.2. 
Exact agreement is shown except for the aerodynamic center location at Mach 1.2 
where the calculated value is 5% low . This is a remarkable agreement and increases 
the confidence in the pressure data for which there are no corresponding wind tunnel 
measurements . 

The magnitude of the flow interference effects is suggested by Table 2-2, which com- 
pares the total normal force coefficient of the mated configuration with that of the 



Figure 2-5. Booster Longitudinal Characteristics in Presence of Orbiter 

2-10 



Table 2-2. Inter-vehicle Flow Interference Effects 
(M = 0.8, a = 0) 


Total Normal Force 
Coefficient 

Booster Alone 

0.0654 

Orbiter Alone 

0.0096 

Total 

0.0750 

Booster in Presence of Orbiter 

0.0925 

Orbiter in Presence of Booster 

0.0083 

Total 

0.1008 


(Interference = 26% of total) 


independent vehicles at 0.8 Mach number. Interference between the vehicles accounts 
for 26% of the total normal force coefficient of the mated configuration. These results 
were obtained from the Woodward analysis. 

Table 2-3 summarizes the total normal force coefficients derived using the Woodward 
analysis at Mach numbers of 0. 8 and 1.2 at 0 and 6 degrees angle -of-attack. The top 
figures correspond to the mated configuration and the numbers in parentheses pertain 
to the booster in presence of the orbiter. The effect of removing the booster aero- 
dynamic surfaces can also be seen. 


Table 2-3. Total Normal Force Coefficient Summary 



M = 

0.8 

M = 

1.2 


With Aero Surfaces 

W/O Aero Surfaces 

With Aero Surfaces 

W/O Aero Surfaces 


0.1008 

-0.0085 

0.0395 

-0.0472 

a = 0° 

(0.0925) 

(0.0072) 

(0.1014) 

(0. 0225) 


0.4921 

0.1342 

0.5359 

0.1753 

a= 6° 

(0.3534) 

(0.0073) 

(0.4211) 

(0.0557) 


( ) = Booster in presence of orbiter. 
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Figure 2-6 shows typical pressure coefficient distributions for the booster body. The 
pressure coefficient is defined by the equation 

ACp = (16) 

T-00 

where p is the local pressure and p m and are the free-stream pressure and dynamic 
pressure, respectively. 

The top portion of the figure shows the effect of booster /orbiter interference. The 
dashed curve represents the pressure distribution for the isolated booster body. The 
corresponding pressure distribution in the presence of the orbiter is shown by the 
solid curve. These distributions are for a line running from nose to tail, 45 degrees 
from the top of the booster. 

The bottom curves of Figure 2-6 show the effect of the booster aerodynamic surfaces 
on the booster body pressure distribution at 135 degrees from the top. All of the in- 
formation shown in this figure is for 6 degrees angle-of-attack and Mach number 1.2. 

Typical wing pressure distributions are shown in Figure 2-7. This figure shows the 
net (upper surface minus lower surface) pressure distributions in the chordwise 



Figure 2-6. Booster Body Pressure Distributions (M = 1.2, a = 6 deg) 
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Figure 2-7. Orbiter Wing Pressure Distribution (M = 0.8, a - 0 deg) 

direction for four span stations. The solid curves are for the orbiter alone and the 
pressures including the presence of the booster are shown by the dashed lines. 

Good pressure data is essential when performing any type of elastic vehicle response 
or stability analysis in which aerodynamic forces are significant. The proceeding 
discussion has presented a brief insight into the complex flow fields associated wito 
vehicles such as space shuttle. The Woodward method appears at this time to be the 
best analytical approach for determining steady-state pressure distributions for com- 
plex configurations at subsonic and supersonic speeds. Further correlations with 
experimental data are required to corroborate these findings. 

2.3 TURBULENCE RESPONSE ANALYSIS 

This section presents a brief review of the equations of motion, a discussion of the 
random and discrete turbulence models used, and details and results of the response 
calculations. A more detailed discussion of the theory is contained in Reference 1. 

Lagrange's equation of motion 


d_ / dL \ _ 9L + Sjfi _ q 
dt\Sqr/ *q r ^q r ^ 
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is used to formulate the response problem where 


L = T - P 
T = kinetic energy 
P = potential energy 
i = dissipation function 
Q r = r^ 1 generalized force 
q r = r ul generalized coordinate 

The deflection (h) at any point on the vehicle is given by 

n 

h(x,y,z,t) =J^h r (x,y,z) q r (t) (18) 

r=l 

where h r is the normalized deflection at the point in the r** 1 mode shape . 

Using normal (with respect to the stiffness matrix) modes of vibration as generalized 
coordinates, the r^ 1 equation of motion is of the form 


M rr4r ^rr^r ^rr ^r ^r 


n 


yi 6 r s M rs ^r Qr 
s=l 


where 


(19) 


M rs = generalized mass 

D rr = generalized damping 

to r = frequency of mode r 

. * I 0,r=s 
rs " ll.r^s 

and the generalized force, Q r , is defined by Equation 4. 

Assuming harmonic motion, Equation 19 can be written in matrix form as 

|q(co)j = - [a(co)] _ 1 J A f ( co)| 


( 20 ) 



where 


A rs = 1 - 6 rs (77) < 1 + i *r> M rs + 

(21) 

_)l,r=s 
rs (o.r^s 

i 

g r '= structural damping coefficient for mode r 
co = gust frequency 

q r = amplitude of the r^ 1 generalized coordinate 


M - Mrs 
iV1 rs 0 

4pb 3 

(22) 

» , . qf Qrf («) 

rf(W> 4 P b3 u 2 

(23) 

Q rs (G0) 

Qrs(W) " 4Pb 3 GO 2 

(23) 


where p is air density, b is reference length, and qj is the turbulence amplitude. 

Transfer functions are formed relating response to turbulence amplitude by the load 
summation method , which is written 

2 n 

H(x,y,z, to) = S F r < x »y> z ) 


1 n _ _ F f (x,y,z,U>) 

+ — > F r (x,y,z, go) q r (U$ + 

q f 


(24) 


where F£ is the response (acceleration, shear, moment, or torque) due to the r** 1 
degree of freedom oscillating at a unit frequency , F r is the load due to the aerodynamic 
pressure distribution developed by a unit amount of the r^ 1 generalized coordinate 
oscillating at frequency go , and F f is the load due to the aerodynamic pressure 
distribution produced by a gust wave of frequency CO. 


For random turbulence, the response power spectral density (PSD), 0 X , is related 
to the gust PSD, 0, and the transfer function magnitude, |H | , by the equation: 


0^60) = JH(CO) | 2 0(60) 


(25) 



The ratio of rms response to rms gust velocity is given by 



1/2 


0 x (CO)dOJ 


( 26 ) 


where to c is the integration cutoff frequency. The "characteristic frequency" of the 
response, or number of zero crossings with positive slope per second, is expressed 
by 



to 2 0x(u>) dco 


1/2 


(27) 


The turbulence PSD applied in this study is the 99^ percentile curve (including the 
dashed portion) shown in Figure 2-8 reproduced from Reference 4. The equation of 
the curve is 


. _ 683.4 (4000.) <28, 

1 + 0.0067 (4000*) 4 - 05 

where E (*) is the turbulence PSD and k is the wave number. As specified by Reference 
4, design turbulence loads are obtained by multiplying the rms loads by a factor of 3. 

In the present study, the spectrum was normalized to produce a unity rms value by 
dividing Equation 28 by its rms value, 1.94 m/s. Design loads and accelerations are 
obtained, therefore, by multiplying their A values by 5.82 m/s. 


Discrete turbulence is simulated by a Fourier series. 

M 

f(t) = a Q + ^ (a m cos co m t + b m sin co m t) (29) 

m=l 

where 


f(t) is the forcing function 
M is the number of terms in the series 
t = time 
UJ m = m tt/T 
T 

a o = f f(t)dt 
-T 
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Figure 2-8. Spectra of Detailed Wind Profiles (Reference 4) 
T 


am =\f f(t) cos co m t dt 


-T 


'm 


4 / 

-T 


f(t) sin co m t dt 


The resulting force time-history f(t) is periodic with period 2T. 
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As prescribed by Equation 29, the force is stationary. It is often more convenient 
to think of the vehicle as being fixed, with the disturbance moving over it. This is 
described by substituting (x/V-t) for t in Equation 29, which yields 

M 

f(t) = a 0 + ]£ a m cos (k m x “ w m l ) + b m sin “ « m t) (30) 

m=l 


where 


V = W m b / V 

x = x/b 

b is the reference length 

V is the vehicle velocity 

x is the distance from a point on the vehicle to the gust reference axis 

By considering that the cosine terms are symmetric and the sine terms are antisym- 
metric , Equation 30 can be written in complex form as 


f(t) = a Q + Re 

M i(w m t-k x) 

2_, a m e 

-Im 

M i(w m t-k x) 

b m e 


m=l 




The vehicle response to this disturbance is given by 


Z(t) = 


M 


a m H R, 


+ b. 


m 


mHl m) 


cos to m t + (b m H 


m 


lR m 


_a m H I 


J sln “m‘] 


(32) 


where H R and Hj m are the real and imaginary components of the m^ 1 transfer function 
evaluated at frequency co m , i.e. , 


H 


m 



+ iHj 

m 


(33) 


The discrete gust used in the present study was the "quasi-square-wave" shape des- 
cribed in Reference 4. The gust amplitude was normalized to 1-m/s. The results 
must therefore be multiplied by a factor of 9 to correspond to the design value specified 
in Reference 4. The gust length was 170 meters, producing a gust period equal to one- 
half the period of the first elastic mode at t = 75 seconds. The Fourier series con- 
tained 40 terms, resulting in the gust shape shown in Figure 2-9. 
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Figure 2-9. Discrete Gust Time History 

A pitch plane stability augmentation system (SAS) was incorporated into the symmetric 
analysis for comparison with the unaugmented vehicle response. The SAS block 
diagram is shown in Figure 2-10. The SAS effects are included by modifying the A- 
matrix in Equation 20. The procedure for doing this is given in detail in Reference 1. 

Propellant sloshing parameters for the booster at t = 60 and 75 seconds are shown in 
Table 2-4. These values were determined by the method of Reference 8. These modes 
were mass coupled to the vehicle modes as described in Reference 1. 



Table 2 

-4. Sloshing Parameters 


Time 

(sec) 

Tank 

Slosh Mass 
(kg) 

Slosh Spring 
(N/m) 

Slosh eg 
(m) 

60 

lo 2 

205,100 

1,162,000 

44.0 


lh 2 

12,720 

72,000 

71.3 

75 

lo 2 

205,100 

1,211,000 

44.8 


lh 2 

12,720 

75,100 

74.4 
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t a 


= 0. 04 sec 


^1 


= 1 

= 0. 35 sec 


Kflg = 0* 35 sec 


SENSOR NO. 1 LOCATED AT BOOSTER NOSE (STATION 25.4m) 

SENSOR NO. 2 LOCATED AT BOOSTER ENGINE GIMBAL (STATION 154.0m) 


Figure 2-10. Pitch Plane SAS 

Responses were computed for the analysis conditions of Table 1-1 for the following 
items: 

1 . Booster wing root shear 

2 . Booster wing root bending moment 

3. Booster wing root torque (about booster station 84.0m) 

4. Booster wing tip acceleration 

5. Booster body acceleration at booster station 29.2m 

6. Booster body acceleration at booster station 44.6m 

7. Booster body acceleration at booster station 59.1m 

8. Booster body acceleration at booster station 76.9m 

9. Booster body acceleration at booster station 99.2m 

10. Orbiter body acceleration at Orbiter station 5.08m 

11. Orbiter body acceleration at orbiter station 24.6m 

12. Orbiter body acceleration at orbiter station 45.2m 

13. Orbiter body acceleration at orbiter station 68.5m 
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Representative transfer functions are shown in'Figures 2-11 through 2-32. All transfer 
functions shown herein pertain to flight at Mach 1.2. (Each figure contains three sets 
of symbols, corresponding to one degree of freedom, total rigid body response, and 
total response including sloshing and elastic modes.) . The one degree of freedom cases 
represent vertical translation and side translation for the symmetric and antisymmetric 
conditions, respectively. 

Figure 2-11 through 2-16 apply to the symmetric flight condition with the SAS off and 
booster aerodynamic surfaces intact (Conditions 1 and 2 in Table 1-1). Corresponding 
plots with the SAS on (Condition 3) are given in Figures 2-17 through 2-22. Antisym- 
metric transfer functions for Conditions 5 and 6 (SAS off, booster aerodynamic surfaces 
on) are contained in Figures 2-23 through 2-28. Transfer functions for the conditions 
having booster aerodynamic surfaces removed are shown in Figures 2-29 and 2-30 for 
tiie symmetric case and 2-31 and 2-32 for the antisymmetric condition. 

The transfer function peaks can generally be seen to coincide with vehicle resonant 
frequencies. In some instances, however, a peak occurs at a frequency well above 
the highest modal value. This can be seen for example in Figure 2-15. These high- 
frequency peaks, probably due to aerodynamic coupling, do not appreciably affect the > 
results , however. In the random case , this is due to the fact that the turbulence 
spectrum drops off rapidly with frequency. Figure 2-33 shows the response power 
spectral density for the transfer function of Figure 2-15. Note that the high-frequency 
peak is now two orders of magnitude below the peak at 2.2 Hz. In fact, the random 
response parameter, X, in all cases had leveled off at 5 Hz. Therefore, 5 Hz was 
used as the cutoff frequency in evaluating A and N 0 by Equations 26 and 27. 

In the discrete analyses, since the gust wavelength was established so as to excite the 
vehicle fundamental structural mode, the high-frequency peak in the transfer functions 
had no appreciable effect on the response time histories . 

Results of the random analysis are summarized in Tables 2-5 through 2-10, Tables 
2-5 and 2-6 give the 7L and N Q values respectively for one degree of freedom. Corre- 
sponding rigid body results are listed in Tables 2-7 and 2-8. The total A and N 0 
values are listed in Tables 2-9 and 2-10, respectively. 

The total A values of booster body and orbiter body acceleration are shown graphically 
in Figures 2-34 and 2-35 respectively for the conditions with booster aerodynamic 
surfaces intact. Figure 2-34 shows that symmetrically the maximum acceleration 
on the booster occurs at 0. 8 Mach number. It also shows the effectiveness of the SAS 
in reducing the accelerations at 1.2 Mach number. The SAS reduces the booster rms 
nose acceleration by 38% and the rms tail acceleration by 7%. In the antisymmetric 
case, the booster body rms side acceleration at the nose is the same at Mach 0.8 and 
Mach 1.2. At the aft end, however, the rms acceleration is higher at Mach 1.2 by 
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Figure 2-11. Booster Wing Root Shear Symmetric 
Transfer Function, M = 1. 2 
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Figure 2-12. Booster Wing Root Bending Moment Symmetric 
Transfer Function, M = 1.2 
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Figure 2-13. Booster Wing Root Torque Symmetric 
Transfer Function, M = 1. 2 
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Figure 2-14. Booster Wing Tip Acceleration Symmetric 
Transfer Function, M = 1. 2 
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Figure 2-15, Booster Body Acceleration Symmetric Transfer 
Function, Booster Station 29. 2m, M = 1. 2 
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Figure 2-16. Orbiter Body Acceleration Symmetric Transfer 
Function, Orbiter Station 5. 08 m, M = 1. 2 
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Figure 2-17. Booster Wing Root Shear Symmetric 
Transfer Function, M = 1.2, SAS on 
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Figure 2-18. Booster Wing Root Bending Moment Symmetric 
Transfer Function, M = 1.2, SAS on 
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Figure 2-19. Booster Wing Root Torque Symmetric 
Transfer Function, M = 1.2, SAS on 
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Figure 2-20. Booster Wing Tip Acceleration Symmetric 
Transfer Function, M = 1.2, SAS on 
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Figure 2-21. Booster Body Acceleration Symmetric Transfer Function, 
Booster Station 29. 2 m, M = 1. 2, SAS on 
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Figure 2-22. Orbiter Body Acceleration Symmetric Transfer Function, 
Orbiter Station 5. 08m, M = 1.2, SAS on 
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Figure 2-23. Booster Wing Root Bending Moment Antisymmetric 
Transfer Function, M = 1.2 
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Figure 2-24. Booster Wing Root Shear Antisymmetric 
Transfer Function, M = 1. 2 
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Figure 2-25. Booster Wing Root Torque Antisymmetric 
Transfer Function, M = 1. 2 
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Figure 2-26. Booster Wing Tip Acceleration Antisymmetric 
Transfer Function, M = 1. 2 
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Figure 2-2 7. Booster Body Acceleration Antisymmetric Transfer 
Function, Booster Station 29.2 m, M = 1. 2 
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Figure 2-28. Orbiter Body Acceleration Antisymmetric Transfer 
Function, Orbiter Station 5. 08m, M = 1.2 
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Figure 2-29. Booster Body Acceleration Symmetric Transfer Function, Booster 
Station 29. 2m, M=1.2, Booster Aerodynamic Surfaces Removed 
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Figure 2-30. Orbiter Body Acceleration Symmetric Transfer Function, Orbiter 
Station 5. 08m, M = 1. 2, Booster Aerodynamic Surfaces Removed 
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Figure 2-31. Booster Body Acceleration Antisymmetric Transfer Function, Booster 
Station 29. 2m, M = 1.2, Booster Aerodynamic Surfaces Removed 
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Figure 2-32. Orbiter Body Acceleration Antisymmetric Transfer Function, Orbiter 
Station 5. 08m, M = 1.2, Booster Aerodynamic Surfaces Removed 
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Table 2-5. 1 Degree of Freedom A Values 



Condition No. 




2 

3 

4 

5 

7 

8 

9 

10 

11 

Booster Wing 

Hoot Shear (10 3 N) 

Root Bending Moment (10^N-m) 
Root Torque* (10^N-m) 

Tip Acceleration (10 -2 g) 

56.84 

57.45 

46.33 

0.937 

56.84 

57.45 

46.33 

0.937 

— 

1.659 

1.861 

5.725 

0.0 

— 

55.66 

61.15 

34.86 

0.483 

— 

2.522 

2.872 

4.325 

0.0 

— 

Booster Body Acceleration (10 _2 g) 
Booster Station 29.2m 
Booster Station 44 . 6m 
Booster Station 59.1m 
Booster Station 76.9m 
Booster Station 99.2m 





O.C 

31 

0.3 

85 

0.0 

21 



0.0 

47 

0.3 

01 

0.0 

54 

Orbiter Body Acceleration (10 _2 'g) 
Orbiter Station 5.08m 
Orbiter Station 24 . 6m 
Orbiter Station 45.2m 
Orbiter Station 68.5m 

0. 

937 

0. 

! 

937 

o.< 

331 

’ 

0.3 

85 

’ 

O.C 

121 

0.' 

f 

183 

O.C 

47 

0.2 

01 

O.C 

64 
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Table 2-6. 1 Degree of Freedom N 0 Values 


Response Item 

Condition No. 


2 

3 

4 

5 

7 

8 

9 

10 

11 

Booster Wing 
Root Shear 

Root Bending Moment 
Root Torque* 

Tip Acceleration 

1.358 

1.467 

1.428 

1.153 

1.358 

1.467 

1.428 

1.153 


2.496 

2.499 

1.326 

0.0 

— 

1,348 

1.390 

1.437 

1.326 

— 

1.763 

1.769 

1.119 

0.0 

— 

Booster Body Acceleration 
Booster Station 29,2m 
Booster Station 44.6m 
Booster Station 59 . lm 
Booster Station 76.9m 
Booster Station 99.2m 





0.2 

85 

0.9 

95 

0.2 

39 



0.2 

31 

0.7 

66 

0.2 

22 

Orbiter Body Acceleration 
Orbiter Station 5.08m 
Orbiter Station 24.6m 
Orbiter Station 45.2m 
Orbiter Station 68.5m 

1 . 

L53 

1J 

L53 

0.2 

85 

0.9 

95 

0.2 

39 

’ 

1.3 

26 

0.2 

31 

0.7 

r 66 

0.2 

22 
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Table 2-7. Rigid Body A Values 


Response Item 


Condition No. 



2 

3 

4 

5 

a 

8 

9 

10 

11 

Booster Wing 










Root Shear (10 2 N) 

61.60 

60.09 

— 

2.881 

— 

55.92 

— 

2.431 

— 

Root Bending Moment (10^N-m) 

64.21 

61.56 

— 

2.431 

— 

61.15 

— 

2.954 

— 

Root Torque* (10^N-m) 

49.30 

49.12 

— 

4.751 

I 

35.80 

— 

4.527 

— 

Tip Acceleration (10 -2 g) 

3.823 

2.480 

— 

1.338 


2.903 

— 

0.363 

— 

Booster Body Acceleration (10“2g) 







■ 



Booster Station 29.2m 

1.480 

0.929 

0.190 

Aiiii 


1.788 

0.098 

0.418 

0.105 

Booster Station 44.6m 

0.698 

0.771 

0.108 


0.026 

0.937 

0.047 

0.367 

0.084 

Booster Station 59.1m 

0.930 

0.994 

0.027 


0.049 

0.465 

0.047 

0.326 

0.063 

Booster Station 76.9m 

2.058 

1.546 

0.083 

0.795 

0.086 

1.323 

0.113 

0.299 

0.036 

Booster Station 99.2m 

3.615 

2.367 

0.213 

1.293 

0.129 

2.715 

0.206 

0.311 

0.012 

Orbiter Body Acceleration (10“ 2 g) 



SSI 







Orbiter Station 5.08m 

3.216 

1.719 

0.337 

0.910 

0.234 

3.385 

0.204 



Orbiter Station 24.6m 

1.867 


0.224 

0.627 

0.195 

2.152 

0.122 



Orbiter Station 45.2m 

0.682 


0.105 

0.578 

0.174 

0.909 

0.045 

SIS 


Orbiter Station 68.5m 

1.499 

1.266 

0.037 

0.979 

0.111 

0.839 

0.079 

0.201 

mil 
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Table 2-8. Rigid Body N 0 Values 


Response Item 

Condition No. 

2 

3 

4 

5 

7 

8 

9 

10 

11 

Booster Wing 


M 




■1 




Root Shear 



— 

B 

— 


— 

1.888 

— 

Root Bending Moment 



— 


— 

Sill 

— 

1.898 

— 

Root Torque* 


1.451 



— 

1.535 


1.123 


Tip Acceleration 

1.330 

1.895 


1.029 

— 

1.550 



1 

Booster Body Acceleration 

mm 



■gt 




■ 


Booster Station 29.2m 


2.476 

0.281 

0.945 

0.336 

2.007 

0.362 

0.496 

0.211 

Booster Station 44.6m 


1.585 

0.280 

0.630 

0.178 

2.172 

0.375 

0.610 

0.215 

Booster Station 59.1m 

1.163 

1.120 

0.291 

0.812 

0.238 

1.375 

0.232 

0.766 

0.221 

Booster Station 76.9m 

1.261 

1.571 

0.296 

0.875 

0.272 

1.301 

0.278 

0.983 

0.243 

Booster Station 99.2m 

1.325 

1.871 

0.289 

0.892 

0.286 

1.536 

0.298 

1.132 

0.331 

Orbiter Body Acceleration 








BB 

m ■ 

Orbiter Station 5.08m 

'mm/ ■ 

2.653 

0.283 

0.779 

0.274 

1.892 

0.344 

0.616 

0.210 

Orbiter Station 24.6m 


2.603 

0.282 

0.920 

0.266 

1.968 

0.356 

0.661 

0.215 

Orbiter Station 45.2m 


1.545 

0.280 

1.075 

0.260 

2.180 

0.374 

0.707 

0.220 

Orbiter Station 68.5m 

1.214 

1.372 

0.317 

1.137 

0.225 

1.052 

0.259 

0.901 

0.243 
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Table 2-9. Total A Values 


Booster Wing 

Root Shear (10 ,? N) 

Root Bonding Moment (l^N-m) 
Root Torque* (lO^N-m) 

Tip Acceleration (1 0 -2 g) 


ltes]x>nse Item 


|Booster Body Acceleration (10" 2 g) 
Booster Station 29.2m 
Booster Station 44.6m 
Booster Station 59.1m 
Booster Station 76.9m 
Booster Station 99.2m 


|0rbiter Body Acceleration (10“ 2 g) 
Orbiter Station 5.08m 
Orbiter Station 24.6m 
Orbiter Station 45.2m 
Orbiter Station 68.5m 


2 

3 

4 

5 

7 

8 

9 

10 

11 

62.45 

56.77 


4.751 


62.06 


3.802 

— 

59.13 

55.28 

— 

5.100 

— 

66.15 

— 

3.666 

— 

50.12 

50.61 

— 

6.334 

— 

46.77 

— 

4.097 

_ 

56.88 

56.77 

— 

4.997 

— 

96.34 

— 

2.080 


8.102 

5.001 

0.227 

2.249 

0.086 

9.609 

0.370 

2.285 

0.303 

3.377 

2.511 

0.145 

1.545 

0.078 

4.626 

0.258 

1.661 

0.229 

5.831 

3.978 

0.080 

1.666 

0.087 

9.080 

0.167 

1.004 

0.158 

4.947 

4.514 

0.101 

1.726 

0.105 

5.813 

0.123 

0.572 

0.054 

6.779 

6.284 

0.217 

1.950 

0.133 

7.551 

0.262 

1.095 

0.096 

28.77 

52.13 

0.521 

3.148 

0.268 

41.86 

0.900 

3.002 

0.399 

4.409 

13.70 

0.275 

1.421 

0.210 

5.388 

0.455 

1.735 

0.275 

4.702 

3.011 

0.146 

1.046 

0.183 

7.587 

0.256 

1.671 

0.214 

13.10 

8.877 

0.093 

1.586 

0.110 

19.05 

0.154 

3.233 

0.037 
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Table 2-10. 


Total N 0 Values 


|Booster Wing 
Root Shear 

Root Bending Moment 
Root Torque* 

Tip Acceleration 


Response Item 


|Booster Body Acceleration 
Booster Station 29.2m 
Booster Station 44 . 6m 
Booster Station 59.1m 
Booster Station 76.9m 
Booster Station 99.2m 


|Orbiter Body Acceleration 
Orbiter Station 5.08m 
Orbiter Station 24.6m 
Orbiter Station 45.2m 
Orbiter Station 68.5m 


Condition No. 


2.006 

2.279 

2.303 

2.550 


2.032 

2.350 

2.235 

2.646 


2.117 

1.786 

2.409 

2.263 

1.905 


1.950 
1.631 
2.000 
2.445 


1.998 

1.633 

2.248 

1.897 

1.645 


1.267 

1.160 

2.132 

2.235 


0.511 

0.490 

0.793 

0.586 

0.413 


0.728 

0.431 

0.608 

1.350 


0.929 

0.885 

1.311 

3.221 


2.041 

1.061 

1.420 

1.493 

1.863 


2.050 

1.351 

0.895 

1.573 


0.465 

0.465 

0.430 

0.376 

0.336 


0.347 

0.298 

0.284 

0.362 


2.086 

1.981 

2.260 

2.279 


2.184 

1.702 

2.621 

2.762 

1.745 


1.939 

1.859 

1.983 

2.717 


0.518 

0.525 

0.777 

0.722 

0.522 


0.720 

0.495 

0.630 

1.309 


10 


1.329 

1.422 

1.151 

1.455 


1.071 

1.014 

0.994 

1.597 

1.274 


1.052 

0.664 

0.948 

1.360 


11 


0.461 

0.459 

0.455 

0.412 

0.476 


0.472 

0.470 

0.471 

0.729 
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Figure 2-34. Booster A Accelerations 
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nearly a factor of 2. Figure 2-35 shows that the orbiter symmetric rms nose 
acceleration is about five times as large as the booster rms nose acceleration. It is 
also interesting to note that while the SAS lowers the booster accelerations, it in- 
creases the orbiter nose acceleration by 82%. 

The booster and orbiter A accelerations are shown in Figures 2-36 and 2-37, res- 
pectively for the conditions with booster aerodynamic surfaces removed. It can be 
seen that removing the booster aerodynamic surfaces lowers the X accelerations by an 
order of magnitude. The booster and orbiter nose accelerations are consistently higher 

at Mach 0.8 than at Mach 1.2. 

Responses due to 1-m/s amplitude quasi-square wave discrete turbulence are shown 
in Figures 2-38 through 2-46. Figure 2-38 shows envelopes of booster body peak 
acceleration for symmetric and antisymmetric turbulence. As in the random case, 
symmetric turbulence produces higher accelerations. Corresponding plots for the 
orbiter are shown in Figure 2-39. Again, the orbiter nose symmetric acceleration 
is considerably higher than for the booster. 

Actual response time histories due to the 1-m/s quasi-square -wave gust are given in 
Figures 2-40 through 2-46. Booster wing root shear, bending moment, and torque 
due to symmetric gust are given in Figures 2-40, 2-41, and 2-42, respectively. 

Booster wing tip acceleration is shown in Figure 2-43. All of these figures show the 
total, one degree of freedom, and rigid body responses separately. Figures 2-44 and 
2-45 show the corresponding booster and orbiter nose acceleration time histories, 
respectively. Total antisymmetric responses are shown in Figure 2-46 for booster 
wing root shear, bending moment, and torque, and booster wing tip acceleration. 

It is interesting to compare the random and discrete turbulence results. Table 2-11 
summarizes the booster wing root shear and bending moment at Mach 1.2 due to 
symmetric random and discrete turbulence. As specified in Reference 4, the rms 
loads from the random analysis based on the 99th percentile PSD were multiplied by 
three to yield design values. The discrete quasi-square -wave gust results were 
multiplied by nine to yield design values corresponding to a 9 m/s amplitude gust. It 
can be seen that the discrete turbulence loads are higher approximately by the ratio 

of the discrete turbulence level to the 3a random value ( 3 ’ indicating that the 

discrete turbulence design criteria may be overly conservative (corresponding to about 
5a random turbulence) . 

The elastic contribution to the total load is also greater for discrete turbulence. This 
is probably because the gust length was set so as to excite the fundamental mode, 
which had a large wing motion. In fact in the random case, the elastic effect was to 
reduce the total wing root bending moment. 
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SYMMETRIC 



Figure 2-35. Orbiter A Accelerations 
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Figure 2-36. Booster A Accelerations (Booster Aerodynamic Surfaces Removed) 
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Figure 2-38. Booster Peak Acceleration Envelopes at 
M = 1.2 Due to 1-m/s Quasi-square Gust 
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ACCELERATION (G) ACCELERATION (G) 
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OR BITER STATION (m) 

Figure 2-39. Orbiter Peak Acceleration Envelopes at 
M = 1.2 Due to 1-m/s Quasi-square Gust 
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Figure 2-44. Booster Body Acceleration at Booster Station 29. 2m 
at M=l. 2 Due to 1-m/s Symmetric Quasi-square Gust 
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Figure 2-45. Orbiter Body Acceleration at Orbiter Station 5. 08m 
at M=1.2 Due to 1-m/s Symmetric Quasi-square Gust 
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Figure 2-46. Booster Total Response to 1-m/ s Antisymmetric 
Quasi-square Gust at M = l. 2 
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Table 2-11. Booster Wing Root Shear and Bending Moment 
Due to Symmetric Turbulence at Mach 1.2 



Shear (10 5 N) 

Bending Moment (lO^N-m) 

Gust Type 

Total 

Rigid 

Elastic 

Total 

Rigid 

Total 

Rigid 

Elastic 

Total 

Rigid 

Random 

3.634 

3.560 

0.074 

1.02 

3.442 

3.737 

-0.295 

0.92 

Discrete 

7.465 

6.416 

1.049 

1.16 

7.316 

6.350 

0.966 

1.15 


To determine the criticality of these loads, booster design wing loads were obtained 
from Reference 9 and are shown in Table 2-12. The maximum aQ loads include a 20% 
allowance for static aeroelastic effects and elastic portion of turbulence response. 

From these comparisons, it appears that including the rigid portion of the gust response 
in the ascent trajectory analysis and multiplying the results by a factor of 1. 2 is con- 
servative for space shuttle, although this should be verified by a static aeroelastic 
analysis . 


Table 2-12. Booster Design Wing Root Shear and 
Bending Moment (Reference 9) 


Design 

Condition 

Total (Rigid + Elastic) 

Elastic 

1 

Bending Moment 
(10 6 N-m) 

Shear 

(10 5 N) 

Bending Moment 
(10 6 N-m) 

Maximum aQ 

28.77 

19.65 

4.77 

3.30 

(Mach 1.2) 





Reentry 

28.94 

18.43 



2 . 5g Maneuver 

23.08 

16.98 



Design Gust 

22.62 

- 

16.53 
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SECTION 3 


CONCLUSIONS AND RECOMMENDATIONS 


A fully reusable space shuttle configuration has been analysed at Mach 
and 12 during ascent flight to determine Its responses to random and discrete atmos 
pheric turbulence. To simulate an expendable booster, the booster aerodynamic 
surfaces were removed for portions of the analysis. 

Aerodynamic interference effects have been shown to be significant for space shuttle. 
The Woodward method appears to be the best analytical abroach for determin^ 
steady-state pressure distributions for complex configurations at sub some and supe 
sonic speeds. The Woodward method has correlated well with wind tunnel to or e 
and moment data for the configuration analyzed. Further testing is desirable Hx> 
provide pressure data for complex configurations such as space shuttle to verify th 

Woodward approach. 

The symmetric random turbulence analysis showed higher accelerations on the booster 
™d 2 at Mach 0.8 ftan at Mach 1.2. Antisymmetrically, booster accelerations 
are higher at 1.2 Mach number. The booster and orbiter accelerations are lower, 
however, by an order of magnitude in the antisymmetric condition than in the sym 

metric case. 

inclusion of the booster stability augmentation system lowered the 

indicating that the discrete gust ^totelw^Xt 1" the discrete 

that elastic response accounted for about 15% of the tota ng gu 

analysis, while in the random case, the elastic effect was small. 

iatic aeroelastie loads, not computed in the present study, are required to determine 
the actual design margin. 

The turbulence response analysis method employed herein is costly and requires 
extensive input data. Simpler methods would be more attractive for 
design and configuration trade studies. The approximate methods would have to be 
evaluated again^the present approach, however, to determine whether the results are 

sufficiently accurate . 
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